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ABSTRACT 

In this paper we present an analysis of the physical and chemical conditions of the 
planetary nebula NGC 40 through spatially-resolved spectroscopic maps. We also in- 
troduce a new algorithm -2d_neb- based on the well-established iraf nebular pack- 
age, which was developed to enable the use of the spectroscopic maps to easily estimate 
the astrophysical quantities of ionised nebulae. The 2d_neb was benchmarked, and 
we clearly show that it works properly, since it compares nicely with the iraf nebular 
software. 

Using this software, we derive the maps of several physical parameters of NGC 40. 
From these maps, we conclude that T e [N n] shows only a slight temperature variation 
from region to region, with its values constrained between ^8,000 K and ^9,500 K. 
Electron densities, on the other hand, have a much more prominent spatial variation, 
as N e [S II] values vary from ^1,000 cm" 3 to ~3,000 cm" 3 . Maps of the chemical 
abundances also show significant variations. From the big picture of our work, we 
strongly suggest that analysis with spatial resolution be mandatory for more complete 
study of the physical and chemical properties of planetary nebulae. 

Key words: atomic data - ISM: abundances - planetary nebulae: individual: NGC 40. 



1 INTRODUCTION 

Stars of low and intermediate mass (0.8 - 8 Mq) rep- 
resent 95% of the evolved stars in the Galaxy and, be- 
cause of the initial mass function, can be important 
sources of enrichment o f He, N and C in the Interstel- 
lar Medium (ISM) (e.g.. lYin et~afll201Ch . Although super- 
novae are expecte d to be the main source of enrichmen t 
for most species ([C hiappini, Ro mano fc Matteuccil |2003), 
Ivan den Hoek fc Groenewegenl (|l997l ) and iMarigd (|200lh 
predicted theoretically the possibility that planetary neb- 
ulae (PNe) progenitor stars also produce O and Ne and 
dredge these elements up to the stellar surface, resulting 
in possible enrichment of the ISM. In any case, PNe and 
the study of their abundances are clearly important tools to 
study stellar evolution as well as -by considering elements 
other than those just discussed- revealing the imprint of 
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ISM abundances when their progenitor stars were formed, 
in the Galactic as well as extragalactic context. 

The outer layers of the intermediate mass stars, which 
were enriched throughout their evolution, are ejected during 
the Asymptotic Giant Branch (AGB) and post-AGB phases, 
and subsequently ionised by the remnant hot nucleus, re- 
ferred to as the central star (CS). Analysis of the ionization 
lines formed in the ejected envelopes of these evolved stars 
is one of the most important ways in which the chemical and 
physical characteristics of the gas are studied. However, the 
great majority of work done in this field utilizes single slit 
spectroscopy to measure the ionization lines and thus is lim- 
ited to obtaining average values across the observed region of 
the slit (or across multiple regions, as it is sometimes done). 
Utilizing single slits is extremely limiting in the analysis of 
these objects, especially considering the high complexity of 
PNe morphology. In single-slit analysis, the possible spatial 
variations of PNe are usually ignored. 

However, spatial variations in the chemical and physi- 
cal conditions of PNe have been detected in previous works 
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using long slits. iGuerrero et al.l {T996), for example, using 
long-slit spectra, found that the N/O ratio varies from 0.4 
to 4.2 from one region to another in the plan etary nebula 
K 4-5 5. Moreover, also using long-slit spectra, iBalick et al.l 
1 19941 ) found N/H abundance gradients in a few PNe. In 
particular they claimed to have found overabundance fac- 
tors of 2 for NGC 6826 and NGC 6543 and 5 for NGC 7009, 
respectively. On the other hand, a nother empirical study 
of NGC 7009 ijGoncalves et al.ll2003h lowered the overabun- 
dance factor to at most 2, and a detailed 3D photoioniza- 
tion modelling of the nebula eliminated completely the need 
of any overabundance to explain the obse rved long-slit as 
well t he narrow-band images of the nebula {Goncalves et al.1 
2006). By inves tigating 13 bipolar PNe usin g long-slit spec- 
trophotometry, IPerinotto fc Corradil {1998) concluded that 
(within the errors) He, O, and N abundances are constant 
throughout the nebulae, and that the Ne, Ar, and S abun- 
dances are also constant, but their face values have system- 
atic increases toward the outer regions of the nebulae. The 
latter was attributed to inaccuracies in the ionization cor- 
rection factors. 

In recent years a new tool for obtaining spectroscopic 
data has been developed. The Integral Field Unit (IFU) is a 
tool that utilizes a square array of lenses remapped to a slit 
and then imag ed through a spectroscope to obtain spatially 
resolv ed data (|Ren fc Allington-Smithll2002l ). iTsamis et all 
(2008) were the first to use this tool to investigate the spa- 
tial variations of PNe parameters. They verified that there 
is variation on the discrepancy of the ++ abundances ob- 
tained from the collisionally excited lines (CEL's) and opti- 
cal recombination lines (ORL's) on the studied PNe. They 
established correlations between the ionization state of the 
gas and the discrepancies of ++ for 2 PNe and between 
the discrepancy of the ++ and the abundances of C ++ and 
++ with the electron temperature for all 3 studied PNe. 

Similar kinds of data can be produ ced by a spec- 
trosco pic mapping technique developed by Mo nteiro et al.l 
(2004), in which multiple parallel long-slit spectra of a 
nebula are interpolated to create emission-line maps. Both 
types of data, those obtained from IFU and those found 
by the spectroscopic mapping technique, enable better 
spatially-resolved analysis of ionised nebulae than those 
based on single long-slit spectra. Spatially-resolved data are 
also critical in constraining photoionization models which 
are now fully three-dimensional and can produce a great 
number of observable quantities, including projected im- 
ages for each emission line. Photoionizati on models that 
utilize these constraint s were obtained by Monteiro et al.l 
l|2004h . iMonteiro etaD l|2005l ). ISchwarz fc Monteiro! (|200rf l 



for NGC 636 9, Menzel 1 and NGC 6781 as well as by 



iMorisset fc Georgievl (2009) 

To further our understanding of the spatial variations 
of the physical and chemical properties of PNe, we present 
in this work results of the study of NGC 40 using spec- 
troscopic mapping . Following the morphological classes of 
IBalick et all ljl987f) . from their Ha, [N n] and [O ill] images, 
the PN NGC 40 has an elliptical shape, which can also be 
seen in the Ha narrow band image obtained in our observing 
run (Figure 1). Two distinct structures are found within the 
main elliptical region: a bright inner rim and an outer shell. 
A filament-like structure can also be identified beyond the 
outer shell, on the northeast region. The [O ill] emission, 



originating in a hotter, low-density medium, is located be- 
tween the hot bubble and the Ha (or [N n]) fil aments, and 
is be tter represented by a ring like structure (|Balick et al.l 

In X-rays NGC 40 also shows a ring-like emission region, 
and the correspondence of this region with the inner rim 
indicates that the bright envelope, seen in the optical image, 
is caused by the shocked CS fast wind. The absence of X-ray 
emission in the region of an apparent apertur e in the inner 
rim s uggests that this wind is quite sph e rical (|Montez et al.l 
120051) . Based on UV data, IClegg et al.1 l| 19831 ) suggest that 
the C IV 1549 A emission from the nebular envelope is too 
strong to have been produced by normal thermal processes, 
and that it is probable that this emission is a consequence 
of processes related to the central star wind. 

We note that there is a discrepancy in temperature de- 
terminations for the CS of NGC 40, a Wolf-Rayet (WC8), 
in the literature. iFeibelmanl (| 19991 ) calls our attention to the 
peculiarity of the fact that NGC 40 is a low-excitation PN 
but seems to be excited by a 90,000 K central star. How- 
ever, a temperature of about 40,000 K, which would be con- 
sistent with a low ex citation PN, is usually adopted for th e 
NGC 40 central star (IStankevicbll983l ; |Pottasch et al.ll2003l ). 
iBianchi fc Grewingj (|l987ft argue that a "carbon curtain", 
located around the CS, could be the mechanism responsible 
for hiding the CS true effective temperature. 

A number of works have dealt with the physical pa- 
rameters and chemis t ry of NGC 40 in the past decade s 
jPottasch et al l 120031 ; iLiu et al.1 l2004al ; iLiu et al.1 l2004bl ). 
These authors investigated different regions of the nebula 
utilizing inconsistent slit configurations or apertures, and 
thus found results that differ greatly from one paper to an- 
other. Even though these works provide a reliable overall 
picture of the object, they fail to give more detailed infor- 
mation that could potentially lead to important clues as to 
shaping mechanisms and enrichment patterns, among oth- 
ers. 

One important practical limitation in dealing with spa- 
tially resolved data is the fact that all currently-available 
tools used to estimate the physical and chemical properties 
of PNe have been limited to single numerical values consis- 
tent with integrated long-slit fluxes. An alternative to the 
above could be to implement scripts written for iraf, which 
would derive the physical quantities in a pixel per pixel, or 
region per region fashion, and then to arrange all the results 
together. However, these kinds of analyses have rarely been 
reported. 

In this work we present a tool, based on the well es- 
tablished iraf nebular package (stsdas. analysis; De Rober- 
tis, Dufour & Hunt 1987), which is capable of working with 
two-dimensional data, such as IFU and spectroscopic maps. 
This tool, the 2d_neb package, has been tested and shown 
to be consistent with results obtained with iraf nebular 
procedures (see Sec. 3.4). We also utilize the 2d_neb pack- 
age to study the PN NGC 40 and present the first spatially 
resolved results of the spectroscopic mapping of this neb- 
ula, for which we have obtained the extinction coefficient, 
electron density and temperature, and ionic abundances of 
several elements. 

This paper is structured as follows. In Section 2 the ac- 
quisition and reduction of the data are discussed. Section 3 
presents the method of data treatment, including a descrip- 
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tion of the spectroscopic mapping technique and explanation 
of those algorithms included in the 2d_neb package that 
were developed to analyse spectroscopic maps. Some bench- 
marks of this new package are also given within Section 3. 
In Section 4 the empirical results for one of the long-slits ob- 
served in this work are compared with those obtained from 
the spectroscopic maps by means of the 2d_neb package. 
Section 5 is dedicated to our spatially resolved results for 
NGC 40, while the summary and conclusions are in Sec- 
tion 6. 



2 DATA ACQUISITION AND REDUCTION 

The observational data used in this work includes both an 
Ha image and optical long-slit spectra. The observations 
were obtained on the night of 2005 October 28 , at the 
2.56 m Nordic Optical Telescope (NOT), located at the Ob- 
servatorio del Roque de los Muchachos (European North- 
ern Observatory, La Palma, Spain). For both the image and 
spectra, the Andalucia Faint Object Spectrograph and Cam- 
era (ALFOSC), which has a pixel scale of 0.189" pixel -1 , 
was used. 

The Ha image has an exposure time of 30 s and the 
seeing during the exposure was 0.89". In this paper this im- 
age is used only for comparison with the Ha spectroscopic 
map that is constructed using the spectra, as described in 
detail below. The reduction of the Ha image included only 
bias subtraction, in addition to flat-field (dome and sky) cor- 
rection. No photometric calibration of this image was per- 
formed. It is worth noting, however, that the resolution of 
the Ha image is better than the resolution achieved with 
the maps. 

For spectroscopy, the long-slit width used was 1.3". The 
ALFOSC was used with the grism number 7, which has 
600 rules mm -1 and a spectral coverage from 3850 A to 
6850 A. The reciprocal dispersion of the binned pixel is 
3.0 A pixel -1 . The slit was located on 16 different paral- 
lel positions across the nebula. The distance between them 
was fixed and set at 3". For 15 of those positions the expo- 
sure times were 3 x 300 s. For slit B only one 300 s exposure 
was taken. Figure [T] shows our NOT Ha image, and also the 
relative slit positions across the nebula. Note, in addition, 
that the slits do not cover the entire object. Slit lengths in 
Figure [T] are reduced for visualization purposes with their 
respective labels used throughout the text. 

The data reduction was performed using the standard 
procedures for long-slit spectroscopy of the Image Reduction 
Analysis Facility (iraf). Bias frames, flat-field, helium- neon 
wavelength calibrations and standard star (G191B2B) expo- 
sures were obtained during the observation runs, and used 
in the reduction/calibration procedure. Three spectra per 
slit position were taken to improve the signal-to-noise ra- 
tio (S/N) and eliminate cosmic rays. The only non-standard 
procedure used was binning on the dispersion axis. This was 
necessary to improve the wavelength (A) calibration due to 
the fact that the helium-neon lamp was resolved in the arc 
frames. 

The final calibrated spectrum for slit G, which passes 
through the central star, is shown in Figure Fluxes here 
are representative of the emission integrated along the slit, 
which we will later call "WN" (whole nebula) when com- 



Ha Image 




Figure 1. NGC 40 image, observed with an Ho filter. The paral- 
lels slit positions are represented by lines and labeled by letters. 
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Figure 2. Spectrum of the NGC 40, integrated over the whole 
slit G. 



paring long-slit measurements and results obtained from the 
spectroscopic maps, in Sections 3.2 and 4. 



3 2D OBSERVATIONS AND TOOLS FOR 
GASEOUS NEBULAE ANALYSIS 

3.1 Spectroscopic Mapping Construction 

Emission-line mapping is a technique in which the spatial 
profiles from a set of parallel slits, for a given emission-line, 
are combined to crea te an emission-line map of the nebula 
ijMonteiro et al.ll2004h . We applied this technique using the 
16 parallel long-slit spectra described above, to obtain the 
emission-line maps of NGC 40. These maps correspond to a 
field of view of 120" x 47" (see Figure 0. 

This method has been previously applie d to Hub- 
ble 5, NGC 6369, Mz 1 a nd NGC 2022 (iRice et all 



20041: iMonteiro et alJ 12004 iMonteiro et~ai1 120051 and 



Mateluna-Perez et al.l [20061 . respectively), providing infor- 
mative quantities such as the total emitted flux of each ob- 
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Ha Spectroscopic Map 
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Figure 3. Ha spectroscopic map of NGC 40, obtained with the 
spectroscopic mapping reconstruction technique. Fluxes are given 
in erg cm -2 s _1 . 



ject and variations of density and t emperature across the 
nebulae, among others. For example. iMonteiro et al.l (|2005l ) 
found indication of the presence of a dense ring and of a 
bipolar structure of lower density for the planetary nebula 
Mz 1. Based on that map, they proposed a tridimensional 
hourglass model to represent the gas distribution of the neb- 
ula. 

In the present paper, the above method was applied to 
generate 31 spectroscopic maps of the PN NGC 40, each cor- 
responding to a different emission line. The generated maps 
were corrected for the effects of field rotation due to the Atl- 
Azimuth nature of the telescope as well a s for differential at - 
mospheric refraction, using the tables of iFilippenkol (|l982T ). 
To avoid the presence of pixels with spurious information, 
all created maps have been submitted to a signal to noise ra- 
tio (SNR) cut. These cuts were done by applying a mask in 
which all the pixels with a SNR lower than a certain value 
were artificially substituted by (zero). For the brighter 
emission-lines, the adopted value was SNR equal to 7. On the 
other hand, this limit did not work for the weaker emission- 
lines, as almost no signal would survive to the noise-mask 
cleaning. To handle this problem, noise-masks based on dif- 
ferent values were tried. Table [TJ presents the list of all the 
31 emission-line maps that were constructed, and the cor- 
respondent SNR cut-offs (SNcut) of these maps. The third, 
fourth and fifth columns of this table show, respectively, the 
integrated spectroscopic map fluxes (F\) -corresponding to 
the emission of the entire observed nebula-, the extinction- 
corrected intensities as will be discussed in Section 5.1), 
and the number of the figure in which the map is presented 
in this paper (Fig.). 

Figure [3] shows the Ha spectroscopic map. As we noted 
in Figure [T] the fact that slits do not cover the entire nebula 
is clearly seen here: the bright edges determine the limits of 
our observations and correspond to the slits A and P. It is 
evident from Figure [3] that the Ha structures of the image 
in Figure [1] are well-reproduced by the spectroscopic map. 
The latter apply not only to the bright elliptical nebular 
components -the inner rim and outer shells-, but also to 
the much fainter structures of smaller scales, as well as to 
the extended filaments. 



Table 1. List of the 31 spectroscopic maps with the cor- 
responding SNR cut-offs. Observed fluxes and intensities 
(extinction-corrected) are normalized to H/3 = 100. 



Line Identification 


SN cut 


Fx 


h 


Fig. 


H10 3797 


7.0 


3.15 


3.98 


- 


H9 3835 


7.0 


6.48 


7.27 


- 


H8+He I 3888 


7.0 


15.40 


17.24 


- 


He+[Ne m]+He I 3968 


7.0 


11.38 


13.79 


- 


[S n] + [S n] 4069+76 


4.0 


0.64 


0.73 


- 


H<5 4101 


7.0 


18.97 


22.43 


- 


H7 4340 


7.0 


38.81 


43.85 


- 


He 1 4388 


2.5 


0.84 


0.05 


- 


He 1 4471 


7.0 


1.75 


1.92 


- 


Mg 1 4563+71 


3.0 


0.15 


0.14 


- 


[C m] + [C iv] 4652 


1.0 


0.14 


0.14 


- 


He 11 4686 


2.0 


0.02 


0.01 


- 


[Ar iv] 4711 


2.5 


0.08 


0.05 


- 


H/3 4861 


7.0 


100.00 


100.00 


5 


He 1 4921 


7.0 


1.03 


1.01 


- 


[O m] 4959 


7.0 


20.38 


19.35 




TO ml ^007 


7.0 


64.22 


58.69 


5 


[N 1] 5198+5200 


5.0 


0.0.45 


0.41 




[CI m] 5517 


4.0 


0.12 


0.10 




[CI in] 5537 


5.0 


0.73 


0.43 




[O 1] 5577 


5.0 


1.01 


0.56 




[N 11] 5755 


7.0 


3.11 


2.57 




He 1 5876 


7.0 


13.61 


11.00 




[O 1] 6300 


7.0 


3.33 


2.56 




[O 1] 6363 


7.0 


0.87 


0.67 




[N 11] 6548 


7.0 


113.54 


84.76 




Ha 6563 


7.0 


400.44 


297.23 


3 


[N 11] 6584 


7.0 


346.77 


253.76 


5 


He 1 6678 


7.0 


3.11 


2.28 




[S 11] 6717 


7.0 


12.41 


9.11 




[S 11] 6731 


7.0 


16.97 


12.41 


5 


Observed H^ flux, F(H/3) 


: 1.52 X 


10 -11 erg 


cm -2 s~ 


1 



3.2 Robustness of the 2D flux measurements 

To check the accuracy of the technique, what follows is the 
comparison between the fluxes measured directly from a sin- 
gle slit and those obtained from the spectroscopic map. Sev- 
eral fluxes from the slit G (see Figure[TJ were compared with 
the fluxes measured from the region corresponding to slit G 
in the maps. 

Flux measurements from slit G were obtained using the 
splot task of iraf. These measurements have been normal- 
ized to correspond to a region whose width is 1 arcsec (fol- 
lowing the fact that the resolution of the maps is 1" pixel - ). 
The comparison is made for 5 separate nebular regions. The 
extent of each of these regions, at the position of slit G, is 
given in Tabled 

The emission-lines used for comparison were selected 
in such a way that isolated, blended, strong and weak lines 
were each chosen. These measurements are shown in Tabled 
together with the relative errors between the two types of 
flux measurements. In more than 75% of the cases, the dis- 
crepancies are lower than 10%. In the worst case, it is 25%. 

The observed discrepancies are likely due to effects spe- 
cific to the regions chosen, such as a single pixel with a dis- 
crepant value, since this value will affect the summed flux 
of the region. Additionally, since the maps are constructed 
by interpolating between the slits, discrepancies might oc- 
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Table 3. Comparison between observed fluxes measured from the slit G and from the corresponding region in the 
emission-lines maps. Fluxes are given in erg cm -2 s — 1 . 



Line Identification 




NOS 






NIR 






SIR 






slit G 


map 


6 (%) 


slit G 


map 


S(%) 


slit G 


map 


8(%) 


H/9 4861 


3.22(-14) 


3.53(-14) 


8.8 


1.18(-13) 


1.39(-13) 


14.5 


1.37(-13) 


1.41(-13) 


2.9 


[O in] 4959 


3.59(-15) 


4.50(-15) 


20.2 


2.80(-14) 


2.86(-14) 


1.9 


3.07(-14) 


3.24(-14) 


5.2 


[O in] 5007 


1.10(-14) 


1.48(-14) 


25.2 


8.49(-14) 


8.85(-14) 


1.0 


9.37(-14) 


9.97(-14) 


6.1 


[N n] 6548 


3.94(-14) 


4.17(-14) 


6.5 


1.36(-13) 


1.60(-13) 


14.9 


1.53(-13) 


1.59(-13) 


3.9 


Ha 6563 


1.30(-13) 


1.39(-13) 


6.2 


4.71(-13) 


5.47(-13) 


13.8 


5.45(-13) 


5.57(-13) 


2.2 


[N n] 6584 


1.18(-13) 


1.25(-13) 


5.6 


4.13(-13) 


4.85(-13) 


14.9 


4.68(-13) 


4.78(-13) 


2.0 


[S n] 6717 


5.10(-15) 


5.24(-15) 


2.7 


1.80(-14) 


1.90(-14) 


5.2 


1.77(-14) 


1.68(-14) 


5.9 


[S n] 6731 


6.62(-15) 


6.87(-15) 


3.6 


2.45(-14) 


2.82(-14) 


13.1 


2.50(-14) 


2.45(-14) 


2.2 



Table 3 — continued 



Line Identification 




SOS 






WN 






slit G 


map 


5(%) 


slit G 


map 


S(%) 


H/3 4861 


3.00(-14) 


3.00(-14) 


0.0 


3.88(-13) 


3.89(-13) 


0.1 


[O m] 4959 


1.98(-15) 


1.84(-15) 


7.5 


7.49(-14) 


7.13(-14) 


5.0 


[O in] 5007 


5.93(-15) 


6.15(-15) 


3.5 


2.25(-13) 


2.21(-13) 


1.8 


[N n] 6548 


3.55(-14) 


3.57(-14) 


0.6 


4.14(-13) 


4.22(-13) 


1.7 


Ha 6563 


1.20(-13) 


1.19(-13) 


0.3 


1.47(-12) 


1.48(-12) 


0.4 


[N n] 6584 


1.08(-13) 


1.08(-13) 


0.5 


1.28(-12) 


1.29(-12) 


0.6 


[S n] 6717 


4.67(-15) 


4.34(-15) 


7.1 


5.25(-14) 


4.64(-14) 


13.1 


[S n] 6731 


5.58(-15) 


5.43(-15) 


2.7 


8.27(-14) 


6.59(-14) 


25.4 



Table 2. Nebular regions of NGC 40, as defined along slit G. 



Region 


ID 


Approximated extension 
(arcsec) 


North Outer Shell 


NOS 


from 46 to 57 


North Inner Rim 


NIR 


from 60 to 76 


South Inner Rim 


SIR 


from 79 to 99 


South Outer Shell 


SOS 


from 101 to 112 


Integrated nebular emission 


WN 


from 46 to 112 



cur due to the fact that the exact slit position is not ex- 
tracted from the map. Figure [4] shows a comparison of the 
spatial profile extracted from the map and that from slit G. 
In this figure the profiles of [O ill] 5007 A, [N n] 6584 A, 
and [S n] 6731 A emission lines are plotted. The profile of 
slit G is represented by the solid line, while that of the spec- 
troscopic map is depicted by the dashed one. Note that, 
with the exception of the central region -where the central 
star is located-, the correspondence between both profiles 
(or emission distribution) is very good. It is also evident 
that the majority of the discrepant regions correspond to 
either high intensity gradients or low SNR regions, in the 
spatial axis. This result shows that the spectroscopic map- 
ping technique is, in fact, a reliable method of representing 
the emission-line fluxes. 



3.3 Tools to deal with 2D emission-line maps 

To perform the usual diagnostic calculations efficiently for 
the emission-line maps, it was necessary to develop a num- 
ber of IDL routines. The set of routines we built compose 
the 2d_neb package, which is described in what follows. The 



software package is divided in two parts: the first one is re- 
sponsible for extinction correction and the second calculates 
the physical and chemical conditions of the nebula (N e , T e , 
ionic and total chemical abundances). The 2d_neb package 
was built based on tasks of the iraf stsdas. analysis nebular 
package, adapted to work with 2D images. As such, all the 
results generated by these tools are 2D maps created in fits 
file format. 

The routines that compose the first part of the package 
are based on the lOsterbrock fc Ferlandl ()2006t ) formalism, 
and are responsible for calculating the extinction coefficient, 
c(H/3), and correcting the data for reddening. The Balmer 
lines are used to determine c(H/3). For these tasks the user 
can choose between Ha, H7 or H5 -relative to H/ 3-, and 
can d ecide which extinction curve will be adopted: ISeatonl 
|l979t ) or Cardelli, Clayton, fc Mathi s (1989), with the coef- 
ficients reviewed bv lO'Donnelll (| 19941 '). When using Cardelli, 
Clayton & Mathis (1989) curve, it is necessary to assume 
a value for the R„ parametrization coefficient, defined as 
R v = (A„ / E(B-V)). The de fault value of R^, u sed for the 
calculations hereafter, is 3.1 (|Cardelli et al.l l 1989V 

The remaining routines of 2d_neb were based on 
the iraf temden task , following the recipes given by 
|Pe Robertis et all (jl987) . These routines are responsible for 
solving the equation of statistical equilibrium. The atomic 
data used in our tools are read from iraf's atomic data 
directory. The electron temperature (T e ) and density (N e ) 
are determined from the solution of the equations of sta- 
tistical equilibrium. This package also calculates chemical 
abundances (ionic and total). The latter procedures were 
based on the ionic and abund tasks of the iraf facility. To 
perform these calculations, the effective recombination coef- 
ficient for the H/3 and the He main transitions, not included 



6 M. L. Leal-Ferreira, D. R. Gongalves, H. Monteiro and J. W. Richards 



[0 III] 5007 A 











■ lorA A , \ 













Table 4. Extinction correction performance. 
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Figure 4. Comparison between the slit and spectroscopic map 
fluxes (in erg cm -2 s — 1 ). From top to bottom it is shown: 
[O III] 5007 A, [N II] 6584 A, and [S II] 6731 A. For each emission, 
the solid line shows the emission measured directly from slit G, 
and the dashed line shows the emission from the spectroscopic 
map. In the bottom of each plot it is shown the ratio between 
both measurements. The region where the difference between the 
fluxes is lower than 10% is highlighted between 2 horizontal lines. 



Line ID (A in A) 


F 


c(H(5) 


l2d_neb 


^redcorr 


[O n] 3727 


33.4 


1.13 


77.2 


77.2 


[O II] 3727 


10.1 


0.70 


17.0 


17.0 


H7 4101 


14.7 


0.26 


16.9 


16.9 


H7 4101 


16.8 


0.43 


21.1 


21.1 


[O in] 5007 


120.8 


0.50 


115.7 


115.7 


[O in] 5007 


1220.0 


0.33 


1185.5 


1185.4 


[N 11] 5755 


2.2 


1.20 


1.3 


1.3 


[N 11] 5755 


3.2 


0.19 


3.0 


3.0 


[S 11] 6731 


20.2 


0.38 


15.3 


15.3 


[S 11] 6731 


17.8 


0.95 


8.8 


8.8 


[O 11] 7330 


48.3 


0.83 


22.5 


22.5 


[O 11] 7330 


30.5 


0.12 


27.3 


27.3 



The fluxes, F, are normalized to H/3 



100. 



in the iraf's atomic data directory, are required. These 
atomic data were taken from Benjamin, Skillman & Smits 
(1999). Total chemical abundance de rivation uses the ionic 
correc tion factors (ICF) formalism of |KiriEsburg h~fc Barlowl 
||1994 ) . 



3.4 2d_neb package: benchmarks 

In this section we show the calculations that were performed 
to benchmark the results obtained from the 2d_neb package. 
This was done by using arbitrary input data covering rea- 
sonable PN values, and the results of each step are checked 
by comparing them with the results given by the tasks of 
iraf's nebular package. 

The first benchmark that we show concerns the perfor- 
mance of the extinction correction routine. As the redcorr 
task (stsdas. analysis. nebular) works with the original coef- 
ficients from the Cardelli, Clayton & Mathis (1989) extinc- 
tion function, we adapted our algorithm to also work with 
these coefficients. We performed the extinction correction 
using this extinction function and adopting R„ = 3.1, as is 
done in iraf. We choose 4 different values for A, well dis- 
tributed inside the optical spectroscopic range. For each of 
these wavelengths, we choose different values of c(H/3) and 
non-corrected observed fluxes. 

Table [4] presents the results of the extinction correction 
for the tested values. The values were used to generate 2D 
images that were then fed through 2d_neb. Table columns 
are, respectively: line identification, observed flux, extinction 
coefficient and the extinction-corrected intensity given by 
the 2d_neb (l2d_neb) and by redcorr (l re dcorr)- Note that the 
last two columns are essentially equal (with discrepancies 
always lower than 1%), indicating that our algorithm for 
deriving extinction corrections is working perfectly. 

The next benchmark we performed concerns N e and T e 
determination. We choose to show the results for typical 
density and temperature diagnostics in photoionized nebu- 
lae, when only the optical spectra is available. For each of 
these quantities we show different values of assumed electron 
temperature or density, and different values for the diagnos- 
tic line ratios. The line ratios were chosen in a way to cover 
the range of possible values given in Osterbrock & Ferland 
2006. T e and N e benchmarking results are shown in Ta- 
ble [5] In this table, columns 1 to 5 show, respectively: the 
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Table 5. T e 


and N e 


benchmarks. 






Output 


Input 


Ratio 


2d_neb 


temden 


8 (%) 




A^e 




T e 


T e 






(cm -3 ) 




(K) 


(K) 




rrt [at 7 1 1 

i e [JN II] 


1,500 


250.00 


7,113.9 


7,114.1 


0.0 


T e [N n] 


2,500 


102.67 


9,397.2 


9,399.2 


0.0 


T e [N n] 


4,000 


70.00 


10 809 8 


10,802.7 


0.1 


T e [0 11] 


1,300 


14.59 


14,287.5 


14,259.7 


0.2 


T 1 tf~^ ttI 

I e [V II) 


3,000 


21.48 


8,271.6 


8,257.8 


0.2 


T e [O II] 


4,300 


24.30 


7,007.5 


7,005.9 


0.0 


T e [O III] 


1,000 


850.00 


7,007.5 


7,008.4 


0.0 


T e [O III] 


3,000 


77.34 


14,222.5 


14,223.0 


0.0 


T e [O in] 


5,000 


301.76 


8,953.0 


8,937.7 


0.2 


T e [S II] 


4,700 


3.73 


9,054.8 


9,056.2 


0.0 


T e [S II] 


2,800 


4.01 


11,878.0 


11,897.0 


0.2 


T e [S II] 


1,300 


5.05 


14,180.2 


14,152.4 


0.2 




T e 






A^e 






(K) 




(cm -3 ) 


(cm -3 ) 




7\r r/~\ ttI 
iV e [0 II] 


7,000 


1.90 


1667.9 


1,667.0 


0.1 


2V e [0 II] 


9,000 


1.48 


1,014.7 


1,014.4 


0.0 


N e [0 II] 


15,000 


2.68 


6,036.2 


6,034.2 


0.0 


JV e [S 11] 


7,500 


0.74 


1,654.9 


1,655.3 


0.0 


JV e [S 11] 


10,000 


0.87 


1,030.8 


1,029.3 


0.1 


iV e [S 11] 


13,500 


0.59 


4,997.2 


4,993.7 


0.1 


2V e [Cl in] 


7,000 


0.91 


3,289.2 


3,290.5 


0.0 


2V e [Cl in] 


12,000 


0.84 


4,777.5 


4,776.7 


0.0 


JV e [Cl in] 


11,000 


1.15 


1,380.8 


1,379.4 


0.1 


Af e [Ar iv] 


9,000 


0.94 


4,952.5 


4,957.8 


0.1 


N c [Ar iv] 


10,500 


1.17 


2,165.0 


2,164.1 


0.0 


N e [Al IV] 


12,000 


1.28 


1,174.4 


1,173.6 


0.1 



quantity (output) to be determined, the assumed T e or N e , 
the line ratio, and the results calculated using 2d_neb as 
well as those obtained by temden (stsdas. analysis. nebular 
— iraf). Inspection of Table [5] clearly shows the similar- 
ity between the two results, which can also be seen in the 
last column, which gives the discrepancies between the re- 
sults from 2d_neb and temden. The discrepancies are, in all 
cases, lower than 1%. 

The last benchmark we show is related to ionic abun- 
dance determination. We choose 4 different ions, and their 
main wavelengths within the optical range, to show the func- 
tionality of our routines. For each ion, we assumed a fixed 
value of electron temperature and density. Table [5] shows the 
results of these benchmarks. The first to the sixth column 
shows, respectively: the ion for which the abundance will 
be derived with its wavelength in A, the adopted electron 
density, the adopted electron temperature, and finally the 
ionic abundance (X/H + ) determined by 2d_neb and by the 
ionic task (stsdas. analysis. nebular). The last colunm shows 
the discrepancies, which are always lower than 1%. 

Some other aspects of the 2d_neb performance were 
also tested, such as the determination of the c(H/3) map, 
the derivation of critical densities, emissivities and wave- 
lengths of each transition, and the calculation of the total 
chemical abundance maps. For space reasons, we not show 
these benchmarks here. 

As a last example of the accuracy of our IDL routines 
in 2d_neb, we make sure that the determination of c(H/3) 



Table 6. Benchmarks of the ionic abundances. Inten- 
sities are normalized to H/3 = 100. 



Ion / A(A) 




2d_neb 
(xlO -7 ) 


ionic 8 (%) 
(xlO -7 ) 


Input: 


N e = 


5,000 cm -3 


T e = 13,000 K 


N+ 5755 
N+ 6548 
N+ 6584 


2.36 
76.01 
235.43 


90.96 
247.7 
260.2 


91.03 0.1 
247.8. 0.0 
260.2 0.0 


Input: 


Ne = 


2,000 cm -3 


Te = 8,500 K 


O 6300 
O 6363 


1.10 
0.70 


40.56 
81.20 


40.54 0.0 
81.16 0.0 


Input: 


We = 


1,500 cm -3 


T e = 13,000 K 


0++ 4959 
0++ 5007 


60.91 
112.77 


273.3 
175.4 


273.0 0.1 
175.2 0.1 


Input: 


Ne = 


5,000 cm -3 


Te = 8,000 K 


S+ 6717 
S+ 6731 


2.76 
9.73 


6.038 
12.12 


6.040 0.0 
12.12 0.0 



agrees with the formulation given in lOsterbrock fc Ferlandl 
(|2006h : 

„™ log[F(\)/F(Hp)) t e - log[F(\)/F(Hf3)] oba 

C{HP) = JixTfiW) ' (1) 

where [F(A)/F(H/?)]t eo is theoretical ratio of another Balmer 
line relative to H/3, [F(A)/F(H/3)] {,,, is the equivalent ob- 
served fluxes ratio, and f(A) is the value of the extinction 
curve in a given A. There is good agreement of the criti- 
cal densities, emissivities and wavelengths, of each transi- 
tion, with the results obtained from iraf's ionic task (sts- 
das. analysis. nebular). The agreement is always equal to or 
better than 99%. 



4 SLIT VERSUS SPECTROSCOPIC MAP 
RESULTS 

In this section we show the results (electron densities, tem- 
peratures and ionic abundances) obtained from a single slit, 
and compare them with those obtained from the correspond- 
ing region of the spectroscopic maps. Following the same 
kind of comparison that we did in Section 3.2, slit G will be 
used. Table [7] lists the extinction-corrected intensities ob- 
tained from slit G, for the six nebular regions under anal- 
ysis: north filament (NF), NOS, NIR, SIR, SOS and WN 
(see Table [2] for the definition of these structures) . Table [8] 
lists the results obtained from these intensities (columns 2 
to 8). The electron densities and temperatures adopted for 
the ionic abundance calculations were N e [S n] and T e [N n]. 
It was not possible to determine T e and N e for the NF. Be- 
cause of this, the N e [S n] and T e [N n] values of the NES were 
adopted to enable the estimation of the ionic abundances of 
the NF region. 

From the results of the slit G physical parameters (Ta- 
ble [8]), one can notice the very low regional variation of 
the T e [N n]. On the other hand, the spatial variation of 
the N e [S n] is significant. These results should be compared 
with those coming from the same region of the spectroscopic 
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maps. The quantities shown with brackets in Table[8] are the 
corresponding mean values from the spectroscopic map for 
the region where the slit G was placed. As we were not able 
to calculate a T e [0 ill] map, the results of this parameter 
refer only to the slit G. 

These results (those obtained directly from the slit and 
those obtained from the maps) are in agreement and, with 
the exception of the N e [S n] estimate of the NIR region (dis- 
crepancy of ~21%), all of them show discrepancy of ~10% 
or less. 

In particular, the spatial N e [S n] variation encountered 
in the analysis of the slit G, can be clearly visualized in the 
N e [S n] map (Figure[7|. The absence of variation of T e [N n], 
from region to region, as shown by the slit G results, can also 
be visualized on the results obtained from the T e [N n] map. 
Note that the histogram of this map is much more concen- 
trated around the mean value than the N e [S n] histogram. 



5 MAPPING RESULTS 

The 31 emission-line maps that were constructed using the 
spectroscopic mapping technique are listed in Table [T] The 
corresponding fluxes and intensities, of the entire nebula 
(WN), as well as the adopted SNR cut-off values (see Sec- 
tion 3.1) are also given in this table. 

The emission-line map of Ha was already shown 
in Figure [3] Other maps -namely: H/3, [O ill] 5007 A, 
[N n] 6584 A and [S n] 6731 A- are shown here, in Figure [5] 



H/3 4861 A 

0.0 5*10"" 1xl0" 16 5xlO" 16 2xl0" ,! 5x10" 15 8xl0" 15 1.1x10"" 1.4x10" 




E , , , i , , , i , , , i , , , i , , , i , , , E 
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orcsec 



[0 III] 5007 A 

0.0 5xl0" 17 IxlO" 16 2xl0"' a 6xl0" 15 IxlO" 15 4xl0" 15 6.5xl0"' 5 9xl0" 15 




20 40 60 80 100 120 

orcsec 



5.1 Extinction Correction 

After creating the emission-line maps, we calculated the 
c(H/3) map using the 2d_neb pack age. Ha, H7 an d H<5 m aps 
(weighted by their fluxes) and the lCardelli et all (| 1 9891 ^ ex- 
tinction curve were all considered in deriving the c(H/3) map. 
The variation of c(H,9) across the nebula can be seen in Fig- 
ure[6] The spectroscopic map clearly shows that c(H/3) is not 
constant throughout NGC 40, and the corresponding his- 
togram shows its dispersion. The mean value computed from 
this map, which includes only valid values (of the pixels that 
survived the noise-mask cleaning) is 0.42. The spatial varia- 
tion of c(H/3), seen in Figure [6] suggests that the amount of 
dust across the nebula is not constant and/or the dus t grains 
have different characterist ics from region to region l|Spitzerl 
Il998l : iMonteiro et al.ll2005l ) . Not e that in studying a portion 
of the nebula IClegg et aTT(|l983l ) found c(H/3)=0.70. In our 
map this region corresponds to an area where c(H/3) values 
are higher than average -with minimum and maximum val- 
ues of ~0.22 and ^1 .03, respectively- whose median, 0.62, 
is i n agreement wit h lClegg et al.l dl983l). 

iPottasch et~ai1 (|2003l ) and lAller et al.l (|l972h . using 
data from various regions of the nebula, found c(H/3)=0.605 
and c(H/3)=0.33, respectively, both in reasonable agreement 
with our map. Studying a low-ionization re gion of the neb- 
ula (not clearly identified in t heir paper). lAller fc Czvzakl 



i| 19791) found c(H/3)=0.65. And lLiu et ail (|2004al ). analysing 
long-slit observations along the NGC 40 major axis found, 
similar to previous results, c(H/3)=0.70. The previous results 
indicate that our mean value of c(H/3) lies in between the 
values found by other authors, who used data from differ- 
ent regions of NGC 40. Furthermore, as we have access to 




3fc 



20 



40 60 
orcsec 



[S II] 6731 A 



0.0 4x10"" 1x10"'° 5x10" 10 IxlO" 13 1.5x10" 15 2x10"'* 3x10"' 




Figure 5. From the top to the bottom, the following emission- 
line maps (in erg cm" 2 s" 1 ) are presented: H/3, [O III] A 5007 A, 
[N ill A 6584 A and fS ill A 6731 A. 
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Table 7. Emission-line intensities of slit G. For WN, whose emission was integrated along the slit, the 
observed flux is also given. Intensities and fluxes are normalized to H/3 = 100. 



Line Identification 


NF 


NS 


Slit Intensities 
NIR SIR 


SS 


WN 
F 


I 


H/3 4861 (crv rm~ 2 s — M 


8 1 3f-1 6") 

O.IOI 1UI 


2 21 f-1 3") 

i.ill lij 1 


8.18(-13) 


9.38(-13) 


2 05C-1 3") 


2 67f-1 2) 




hi n ^7Q7 

11 1U Ol £7 1 




4.80 


4.84 


5.68 


5.56 


5.34 


6.84 


He i 3820 






0.85 


1.86 




3.11 


3.97 


HQ 

11C7 JOiJiJ 




7.45 


7.59 


8.87 


7.18 


8.49 


10.80 


fNe tttI 3869 

1 1 1 t, 111J iJOUC/ 






1.05 


- 




1.19 


1.50 


H8+Hp t 3888 

llu 1 ll*' 1 yUUO 




15.51 


17.85 


21.57 


15.8 


22.94 


28.90 


Hf-t-fNe iiiI-I-Hp i 3968 

1 lf_ |^ 1 ^ t. 111 1^ ± It. 1 iJt7UO 




15.88 


15.39 


16.37 


15.71 


16.44 


20.79 


TFp ml +Hp t 4flOQ 

It, HI 1^ lit, 1 rrtJt/f 






0.18 


- 








Hp t 4D9fi 

lit; 1 '-tVJ^XJ 






1.49 


1.17 




1.39 


1.70 


S it 4033 






0.33 


- 








[c; rrl-l-f^ ttI 4f)fiQ-J-7fi 

[O 1 1J [O 1 1J 4:UUjt ' U 




2.24 


2.74 


4.48 


2.46 


9.91 


12.01 


H<5 4101 




25.28 


24.70 


25.54 


24.41 


19.34 


23.27 


H7 4340 




44.85 


45.55 


46.96 


45.42 


40.99 


46.51 


rn T tt1 zisfi'H 






0.70 


0.90 




2.66 


3.00 


He 1 4388 






0.36 


0.53 








He 1 4438+? 






- 


1.09 




9.49 


10.49 


He 1 4471 




2.67 


2.90 


3.18 


2.39 


3.89 


4.27 


N m 4524 

IN 111 rrt^^rt 






- 


0.18 








He tt 4542 

lit- 11 t -±.\J'~T.Zj 






- 


0.32 








Mo- 1 4563-1-71 




1.01 


0.51 


0.58 


0.75 






TC! TTTl + fC ivl 4652 






2.38 


12.27 








Hp tt 4686* 

11L 11 rtUOU 






0.70 


3.88 








\ Ar ivl 471 1 

Al 1 V rrl 11 






0.51 


0.56 


0.78 






H/3 4861 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


100.00 


Hp t 4921 

11C 1 rt£/i*l 






1.20 


1.43 


0.89 


1.71 


1.69 


rn tttI 4Q5Q 




10.91 


23.01 


21.98 


6.49 


18.76 


18.38 


rn tttI ^007 

[W 111J JUU I 




33.14 


68.98 


66.25 


19.2 


57.1 


55.39 


111 II JlfO |^ i. J V J t/ 




1.37 


0.92 


0.93 


1.34 


0.33 


0.31 


[X t_ 11J -j- l_ 11 1J JZU^T ' u 






0.41 


0.46 




\ 4 


1.29 


TPl tttI ^1 7 

111J JJ1 1 




0.26 


0.34 


0.37 


0.44 


0.58 


0.52 


TPl tttI ^^^7 

[v^l 111J OOO 1 




0.28 


0.37 


0.32 


0.47 


0.56 


0.50 


[O 1] 5577 


- 


0.90 


0.41 


0.37 








[N n] 5755 




2.97 


2.92 


2.99 


2.92 


3.51 


3.02 


He 1 5876 




9.60 


11.61 


12.38 


9.11 


20.62 


17.51 


? 5891 






0.47 


0.47 




4.23 


3.59 


[O 1] 6300 




5.47 


3.73 


3.48 


5.11 


4.78 


3.87 


[S m]+He II 6312 




0.54 


0.74 


0.69 


0.23 


1.03 


0.83 


[O 1] 6363 




1.78 


1.22 


1.19 


1.63 


1.42 


1.14 


He II 6406 






0.12 






1.16 


0.93 


[N 11] 6548 


110.20 


87.57 


84.10 


85.19 


88.53 


105.90 


83.45 


Ho 6563 


310.62 


296.92 


296.43 


300.18 


298.7 


377.43 


296.95 


[N n] 6584 


375.66 


267.88 


258.05 


256.87 


269.22 


328.89 


258.18 


He 1 6678 




3.04 


3.24 


3.59 


2.96 


4.00 


3.11 


[S n] 6717 


34.38 


11.60 


11.28 


9.57 


11.45 


13.09 


10.13 


[S n] 6731 


18.42 


15.03 


15.34 


13.47 


13.66 


20.77 


16.05 



* The He II 4686 emission-line of SIR is probably blended with O II 4676 emission-line. 
We were not able to give a clear identification to the lines marked with "?" . 



the spatial variation of c(H/3), it is clear why previous works 
found different extinction constants. The trend of deriving 
higher c(H/3) values in particular portions of the nebula is 
clearly identified in the map of Figure [B] 

The c(H/3) map was used for extinction correction of 
the 31 emission-line maps. In this process, each pixel of each 
spectroscopic map had its extinction correction done based 
on the value of its corresponding pixel in the c(H/3) map. 
This procedure was performed by 2d_neb. 

After the extinction correction we calculated the mean 
intensities for each emission-line of NGC 40 by integrating 



the emission of each spectroscopic map. These mean intensi- 
ties were normalized to H/3 = 100 and are shown in the last 
column of Table [1] Section 3.1. Note that, as a consequence 
of using the c(H/3) map to correct each emission-line map, if 
a given line-map pixel has survived the noise-mask cleaning 
in its map, but has not survived this process in the c(H/3) 
map, then the pixel is not a valid pixel in the subsequent 
extinction-corrected emission-line map. As such the mean 
intensities, shown in the last column of Table [1] may not 
correspond to the same number of pixels as the mean fluxes 
(shown in the third column of the same table). 
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Table 8. Mean electron densities, temperatures and ionic abundances: slit G versus spectroscopic mapping results. 



Parameter 


NF 


NS 


Slit G results 
NIR SIR 


SS 


WN 


Map results 




Literature 
















Mean 


Inf.L. 


Pottasch 1 


Liu 2 


N e [S II] (cm" 3 ) 
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- 
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- 
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9.21(-7) 
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7.74(- 
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C1++/H+ 5517 
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1.88(- 
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Pott asch 1 : iPottasch et alj ll2003h . 
Liu 2 : iLiu et alj l|2004al lrj). 

*The central star was not considered in the determination of the physical and chemical parameters of whole nebula, WN. 



5.2 N e and T e Physical Properties 

Using the 2d_neb the T e [N n] and N e [S n] maps were cal- 
culated. For N e [S ii], we first assumed a constant elec- 
tron temperature, equal to 9,000 K, throughout the map. 
The resultant map was subsequently used to calculate the 
T e [N n] map, and then the N e [S n] was recalculated using 
this T e [N ii] map. 

The region with valid pixels in the T e [N n] map was lim- 
ited by the noise-mask cleaning of the [N n] 5755 A emission- 



line map. As the determination of N e [S n] depends on the 
electron temperature, we assumed a constant temperature 
equal to 9,000 K for the non- valid pixels of the T e [N n] map. 
T e [N ii] and N e [S n] maps are shown in Figure [7] In the 
N e [S ii] map, the dashed line indicates the outer limit be- 
tween the region in which N e [S n] has been calculated by 
adopting T e [N n] = 9,000 K, and the region in which the 
T e [N ii] map was used. 

The mean values of T e [N n] and N e [S n] are, respec- 
tively, 8,850 K and 1,650 cm -3 . Comparing with the liter- 
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Figure 6. Top: c(H/3) map of NGC 40. Bottom: The histogram 
from the c(H/3) map. The mean value of the c(H/3) map is 0.42. 
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ature (Table [S] columns 10 and 11), the mean values we 
found for N e [S n] and T e [N n] are marginally in agree- 
ment with previous determinations of the physical proper- 
tie s in NGC 40. Notice that the N e [S n] value determined 
bv lPottasch et alj l|2003h (2,100 cm" 3 ) is very similar to the 
values that we found for the NIR and SIR regions, when we 
used the data from the slit G (and from the N e [S n] map, 
at that same region). On the other hand, one can notice 
that the N e [S n] map shows great variations, even if we con- 
sider th e internal region only. The values found bvlLiu et alj 
|2004al ) (1,750 cm -3 ) and bv lDelgado-Inglada et all l|200g| ) 
(1,800 cm" 3 ) approach our mean N e [S n] map result. They 
both used the same line intensities . The T e [N n] results 
fromlLiu et al l l|2004al ) (8,400 K) and lDelgado-Inglada et ail 
l|2009l ) (8 ,600 K) again are closer to the mean value of 
the T e [N n] map than the one found by iPottasch et al.l 
(2003). Actually, when we look at the T e [N n] map, we 
hardly see any region with values as low as that found by 
IPottasch et all (|2003l ) (7,500 K), and most of the nebula 
shows T e [N II] > 8,000 K. 

It was not possible to calculate the T e [0 in] spectro- 
scopic map due to the very faint emission of [O ill] 4363 A. 



Figure 7. From the top to the bottom, the following NGC 40 
physical properties plots are shown: T e [N n] map, (in K), 
N e [S II] map (in cm" 3 ), T e [N n] histogram and N e [S II] histogram. 
The region inside the dashed line, on the N e [S n] map, was calcu- 
lated adopting the electronic temperature from the T e [N n] map. 
A temperature equal of 9,000 K has been adopted for the outer 
region. 



In cases like this, the usual analysis made using a single slit 
has an advantage since the emission is integrated through- 
out a region, or throughout the whole nebula, thus the very 
faint emission-line flux is improved in signal to noise. This 
procedure may allow the calculation T e [0 in]. In fact, in 
the case of NGC 40, we did derive T e [0 lll]in two nebular 
regions (NIR and SIR) using the fluxes of slit G, as shown 
in Table U 



5.3 Ionic and Total Chemical Abundances 

T e [N n] and N e [S n] maps were used to calculate the 
ionic abundance maps. For the non-valid pixels of these 
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maps we assumed constant values, T e = 9,000 K and 
N e = 1,600 cm -3 . When more than one emission-line was 
available for determining the ionic abundance of a certain 
ion, all the strong lines were used, each resulting in an inde- 
pendent ionic abundance map. Then, the ionic abundance 
map was generated from these individual maps weighted by 
intensity. These are the ionic maps shown in Figure [8] and 
Figure [9] 

The mean values of the ionic abundances are listed in 
Table [5] in the first column of the map results. We draw 
attention to the fact that these mean values were calcu- 
lated only over the pixels that have survived the noise-mask 
cleaning. This means that the abundance maps which have 
non-valid pixels inside the nebular regions whould have, in 
reality, lower mean abundances if the entire nebula were 
considered. 

As an attempt to determine the abundances represen- 
tative of the whole nebula (WN, in the tables) for the cases 
that have non- valid pixels inside the nebular region, we cal- 
culated, for each abundance map, the sum of the ionic abun- 
dances of all pixels and divided the result by the number 
of pixels that compose the nebulaQ. These results represent 
the lower limits of the corresponding NGC 40 abundances. 
These lower limits (Inf.L.) are also shown in Table [8] col- 
umn 9. 

The mean ionic abundances determined in this paper 
can be compared with the results found in the literature, by 
inspecting Table [8] For some cases the lower limit (column 
9) is closer to the previous ionic determinations. That is 
obvious for He ++ /H + A4686 abundance, whose map shows 
almost no presence of that specie. The ionic abundances 
from C1++/H+ A5517, C1++/H+ A5537, N+/H+ A6583, 
0++/H+ A4959, 0++/H+ A5007, S+/H+ A6717 and 
S + /H + A6731 show a somewhat different result when com- 
pared with the literature. The discrepancies vary from 70% 
up to a factor of 3. These discrepancies can be due to 
the fact that here, we are analysing a bigger area of the 
object: the whole nebula, in stead of one long slit, as in 
the case oflLiu et al. 1 (l2004ah . and a set of small portions, 
in the case of iPottasch et alj (|2003|)). Notice, for instance, 
that despite the fact that the mean value of our Cl ++ /H + 
map (3.9 9 x 10 ~ 7 ) is higher than the mean valu e from 
iLiu et al.1 |2004al ) (4.19 x 10~ 8 ) and lPottasch et al.1 (120031 :) 
(5.9 x 10 -8 and 7.5 x 10~ 8 ), the inner rim presents a lower 
Cl ++ /H + abundance when compared to the shell (see Fig- 
ure [HJ. The result we found for the inner rim is closer to the 
results from the literature. 

Examining the ionic distributions of 0°/ + H+, N+/H+, 
S+/H+, He+/H+ and that of 0++/H+, in Figure [8] and 
Figure |9j we see a region (a diagonal, from left to right, 
that includes the location of the central star) in which ionic 
fractions coming from the low-ionization species (0 + /H + , 
N + /H + and S + /H + ) are less abundant when contrasted 
with the neighbor-nebular portions. On the other hand, the 
map of ++ /H + and He + /H + show values higher than the 
average at the same region. It is worth noticing that, even 
though He + /H + is a low-ionization specie, a great amount of 



1 Since all abundance maps are restricted by the locations where 
the c(H/3) map is valid, we used the number of valid pixels of the 
c(H/9) map for the number of pixels that compose the nebula 
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Figure 8. From the top to the bottom, the following ionic abun- 
dance maps are presented: Ar IV, CI III, He II and O I. 
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He is also expected to be present in this low ionization neb- 
ula. In spite of the peculiar behaviour of the low-ionization 
ions, as well as ++ /H + and He + /H + , the electron temper- 
ature maps (particularly relevant for abundance determina- 
tions) have no obvious gradient (see corresponding map in 
Figure that would justify this result. 

Is this higher ionization tunnel of emission in NGC 40 
similar to the jet structures of a number of other nebulae? 
In NGC 7009, for instance, two pairs of low-ionization knots 
are clearly seen, in addition to a pair of jets, in which the 
excitation degree is significantly higher than that of the im- 
mediate vicinity (|Goncalves et al]|2003l ). However, contrary 
to what we find for NGC 40, in the case of NGC 7009 the 
electron density of the higher-ionization emission tunnel is a 
factor of 2 lower than in the other nebular regions. But note 
that in both cases the T e is roughly the same throughout 
the nebula. 

Concerning the element abundances of helium in 
NGC 40, He/H, we found the mean value of 7.08 x 10" 2 . 
This result has been obtained by summing the He + and 
He ++ ionic abundance maps, and taking its mean value. 
This extremely low total helium abundance can be explained 
by the fact that a significant amount of He is expected to 
be present in a low excitation nebula like NGC 40, and He 
was not considered i n the He/H given ab ove. To circumvent 
this kind of problem, IZhang fc Liul (|2003h suggested that an 
ICF based on the S + and S ++ abundances can be used to 
account for the He abundance of low excitation PNe. Using 
our own data we were not able to calculate a S ++ abundance 
m ap, nevertheless, if we c onsid er the S ++ abund ances found 
by IPottasch et all (|2003l ) and iLiu et ail (|2004bl ) (Table EJ), 
and the mean value of the S + /H + map that we calculated 
(8.88 x 10 -7 , in Table|5}, we can determine a more reliable 
value for the He total abundance of NGC 40. Doing this exer- 
cise, we find He/H {/CF _ P) =9.32 x 10~ 2 , adopting S++/H+ 
from Pottasch et al. (2003), and He/H (JCF _ i) = 1.18 x 10" 1 
by using S++/H+ from Liu et al. (2004b). These three He/H 
are shown in Tabic [8] Notice that t he int erval of S ++ abun - 
dances from IPottasch et al.l i^OOSt ) and ILiu et al.l 11200401 ) 
contains the mean S ++ abundance we found from our slit G, 
if we consider the NS, NIR, SIR and SS regions (again, that 
is not the case for the WN region, because of the contamina- 
tion from the central star region). And, finally, not only we 
believe that 9.32 x 10~ 2 and 1.18 x 10 _1 should correspond 
to the lower and upper He abundance limits of NGC 40, 
but they are also similar to the previous results found by 
ILiu et all l|2004bl ), who adopted the same ICF scheme to 
account for the unobserved amount of He in the total he- 
lium abundance of NGC 40. 

The total oxygen abundance map is shown in Fig- 
ure llOl This map has been obtained summing the ionic abun- 
dance maps of O and 0++ -since, unfortunately, + map 
is not available. This element abundance is, in fact, only 
a lower limit. The corresponding mean value of the O/H 
map is 4.10 x 10~ 5 ; a much lower value when compared 
with the results from the literature. That is a consequence 
of the absence of th e Q + abundance map . As shown by 
IPottasch et all (|2003h and lLiu et ail (l2004bh . this ion dom- 
inates the oxygen abundance of this nebula. 



Figure 9. From the top to the bottom, the following ionic abun- 
dance maps are presented: He I, N II, O III and S II. 
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Figure 10. Inferior limit of the oxygen chemical abundance map. 



6 SUMMARY AND CONCLUSIONS 

In this paper we presented a spatially-resolved analysis of 
the physical and chemical characteristics of the planetary 
nebula NGC 40. To reach these results, 16 long-slit parallel 
optical spectra were obtained, from which we created 31 
emission-line spectroscopic maps. These maps provide the 
spatially-resolved characteristics of the line emission, and 
the total fluxes for the whole nebula. 

Using the usual diagnostic ratios, these maps were ap- 
plied to calculate the 2D density, temperature and abun- 
dance information for this nebula. To derive these quantities, 
we adapted t he tasks from the well-e stablished iraf neb- 
ular package (|De Robertis et ai]|l987t ). to work with two- 
dimensional data, thereby developing, and proving the effi- 
ciency of, the 2d_neb package. 

Each emission-line map was corrected for reddening, 
pixel by pixel, using the c(H/3) map that was created from 
the Hq/H/3, H7/H/J and H5/H/3 map ratios. The c(H/3) map 
shows interesting features, with values that vary from less 
than ~0.2 up to more than ~0.6, as shown in the histogram 
of Figure [6] with a mean characteristic of the whole nebula, 
c(H/3) equal to 0.42. The above spatial extinction variation 
could indicate that dust is mixed with the ionized gas of 
NGC 40. 

Our T e [N 11] map shows only a slight variation from 
region to region, from ~8,000 K up to ~9,500 K. On 
the other hand, these variations are much less prominent 
than those found in the N e [S n] map, where the values 
varies between ~1,000 cm -3 and ~3,000 cm -3 . The his- 
tograms of both physical parameters, T e [N n] and N e [S n] , 
show this dispersion. The density map also shows that 
there is a trend to lower density regions on the west 
side of the nebula, but it does not show any trend for 
the outer regions compared to the inner ones. The mean 
N e [S n] (1,650 cm" 3 ) and T e [N n] (8,850 K), from the spec- 
troscopic maps, are in good agreement wi th other results 
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from the literature dDelgadoTnglada et al.ll2009 | : 
2004al : IPottasch et all I2003L IClegg et all Il983l : (Aller et al 
19723 ). On the other hand, the maps show spatial variations 
that were not addressed by previous works. 

Several of the ionic abundance maps show significant 
spatial variations, i) The Ar 3+ /H + and Cl ++ /H + maps show 



a ring structure on the inner region of NGC 40, and also high 
concentrations from these species on the outer region of the 
nebula, ii) The O /H + e ++ /H + maps show strong spa- 
tial variations. In the first case, the contrast of the regions 
with O /H + < 10~ 5 is clear. For the second, most of the 
nebula has ++ /H + between 10" 5 and 3xl0~ 5 . However, 
regions with lower ++ abundances are present in the outer 
regions of the nebula, and regions in which ++ /H + varies 
from ~5xl0~ 5 up to ~1.3xl0 -4 are found, primarily be- 
tween the inner rim and outer shell, iii) The He + /H + abun- 
dance is also inhomogeneous. iv) The N+/H+ abundance 
map shows that almost the whole nebula has an N + /H + 
between 5xl0 -5 and 9xl0~ 5 or that N + /H + does not vary 
significantly throughout the nebula, v) Finally, the S + /H 
map shows clumps, in which the S + /H varies between less 
than 4xl0~ 7 up to more than 10~ 6 . 

It is important to point out that maps represent a pro- 
jection of the nebula along the line of sight on the sky, and 
so nothing can be said about the tridimensional state of the 
gas (T e , N e , ionic as well as total abundances) from our 
analysis of NGC40. 

However, the results obtained in this work allow us to 
conclude that the spatial variations found in the c(H/3), 
T e [N 11], N e [S 11], and abundance maps unequivocally con- 
firm the necessity of an analysis with spatial resolution for 
a more complete study of the physical and chemical prop- 
erties of planetary nebulae. Furthermore, the spectroscopic 
mapping technique used in this paper and in previous works 
shows important advantages when compared with the usual 
technique that uses one or a few long-slits to analyse the 
obje ct. As an extra ingr e dient in this discussi o n, we note 
thatlCorradi et all (Il997t). IPerinotto fc Corradl (|l998h . and 
iGoncalves et all (|2003l ,"l2006) did not find significant spatial 
variation of the chemical abundances when analysing a num- 
ber of other PNe. As we do find these variations in NGC 40, 
an important question that remains to be addressed is: if 
analysed with the spectroscopic mapping technique, as in 
this work, would those ~20 PNe studied with single long-slit 
spatially resolved spectra (above references) show significant 
chemical variations? 

Also in the context of distant extragalactic PNe, it is 
clear that spatially resolved data using techniques similar 
to the one applied here to study NGC 40 are desirable. 
However, our telescopes/instruments prevent us of such a 
detailed study of distant objects. Fortunately in these situ- 
ations most of the emission of a PN is captured in a single 
slit observation -making the average values more representa- 
tive. And, therefore, the interpretation of results for distant 
extragalactic PNe are, at least in a statistical sense, not as 
seriously affected by the lack of spatial resolution. 
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